Remarks on the non-equilibrium effects and collision dynamics in heavy-ion collisions 
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We study the beam energy dependence of equilibration process and space-time characteristics of 
participant and spectator matter. For this, we simulated the semi-central collisions of 40 Ca + 40 Ca 
at incident energies of 400, 600 and 1000 AMeV within the quantum molecular dynamics (QMD) 
approach. Our numerical calculations based on the molecular dynamics approach show that incident 
energy of the projectile influences the reaction observables drastically. The effect is more visible for 
transverse expansion of the nuclear matter and transparency behavior. The degree of thermalization 
of participant matter, however, remains independent of the incident energy. The characteristics of 
the trajectories followed by the nucleons suffering maximal and minimal binary collisions are also 
analyzed. 

PACS numbers: 25.70.-z, 25.75.Ld, 24.10.Lx 



I. INTRODUCTION 

The simulation of heavy-ion collisions at intermediate 
energies has always played apivotal role in understanding 
the reaction mechanism [l|-|7[ , nature of nucleon-nucleon 
(n-n) interactions, as well as thermalization achieved 
by the nuclear system jsl— f3~ol) - Similar efforts are also 
made at lower tail of the incident energy where fusion 
and cluster-decay processes are dominant (TT| - [T3 |. The 
choice of various transport codes, beyond doubt, pro- 
vides a unique platform to study the properties of hot 
& dense nuclear matter formed during these collisions. 
In these transport models, one starts from the well sep- 
arated projectile and target nuclei and then evolves to- 
wards non-equilibrium stage due to the overlapping of 
two Fermi spheres with large relative momentum. The 
dynamics of heavy-ion (HI) collision can be interpreted 
in terms of useful reaction observables such as collective 
flow [l5l - [l8| . multifragment-emission [3, HE Esj] , tempera- 
ture and entropy production pj} , momentum anisotropy, 
subthreshold particle production as well as elliptical flow 

seeeeiii. 

The different studies have shown that for a given 
colliding geometry, the final particle abundances are 
sensitive towards the incident energy of the projectile 
[23l [24|]. These theoretical predictions, however, were 
largely based on the assumption of thermal equilibrium 
of the nuclear system. It may be emphasized that at SIS 
energies, particles are produced in highly non-equilibrium 
situation much different from the saturation characteris- 
tics of the normal nuclear matter [25| . Recently, the BUU 
simulation of 197 Au+ 197 Au collisions at 1 AGeV reported 
that spectator matter evolution is strongly influenced by 
the compression and expansion of the participant mat- 
ter [26l | . The transverse expansion of the spectator mat- 
ter is found to be sensitive towards the nuclear incom- 
pressibility of the participant zone [26j |. At relativistic 
bombarding energies, the spectator matter fragmentation 
[E Q and its universality characteristics 0, [27l - l29| have 
been extensively studied in recent literature. As far as 
particpant matter physics is concerned, extensive study 



has been made in the past to infer the properties of hot 
and dense participant zone @, H, H^-Hlj. Unfortunately, 
excited nuclear matter under the extreme conditions of 
temperature and density exists for a very short duration. 
This makes it difficult to extract information about the 
nuclear forces and reaction dynamics in experiments. For 
instance, nuclear matter at nuclear density p >1.5p lasts 



34]. It is inter- 



for the time span less than 20 fm/c [22 
esting to study the behavior of participant matter and 
its response towards the evolution of the spectator mat- 
ter during the collision process. This study also bears 
relevance due to its relation with the nuclear equation of 
state (EoS) and flow systematics. 

During the collision process, different kinds of inter- 
actions are at work which are important in their own 
|35l l36j . One is also interested to know behavior of nucle- 
ons in the presence of other nucleons in the surroundings. 
From the trajectories traversed by the nucleons in phase 
space, one can gain important information on the nature 
of the hadronic forces as well. In the present work, we 
aim to analyze the reaction dynamics in the semicentral 
collisions of 40 Ca + 40 Ca at incident energies of 400, 600 
and 1000 AMeV respectively. We shall also examine the 
behavior of nucleons facing maximal and minimal num- 
ber of collisions via stopping and trajectories followed in 
the phase space during the collision process. 

For the simulation of heavy-ion reactions, we shall 
utilize the qua ntum molecular dynamics (QMD) model 
[E IE IE EE HEUll as primary transport theory which is 
described in |TTJ Section IIIII deals with the calculations 
and illustrative results, which are summarized in section 

Uvl 
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THE MODEL 



The quantum molecular dynamics is an N-body the- 
ory that simulates the heavy-ion reactions at intermedi- 
ate energies on event by event basis. This is based on a 
molecular dynamics picture where nucleons interact via 
two and three-body interactions. The explicit two and 
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three-body interactions preserve the fluctuations and cor- 
relations which are important for A-body phenomenon. 
In this approach, each single nucleon of the two colliding 
nuclei is described by a Gaussian wave packet in mo- 
mentum and coordinate space. The centroids of these 
Gaussian wave packets propagate in coordinate (7^3) and 
momentum (V3) spaces in accordance with the classical 
equations of motion: 

dUi . Pi dUi ,„ . 

Pi = 5I~) r i = / „ . + -KT- I 1 ) 
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In the quantum molecular dynamics approach, nucleons 
move along classical trajectories and undergo stochas- 
tic scattering. Stochastic means that scattering ampli- 
tude doesn't relate the scattering angle with the impact 
parameter in a unique way like the collision of classical 
billiard balls. Collisions are Pauli blocked if scattered 
nucleons enter the phase space region which is already 
occupied. The local Skyrme interaction used here, con- 
sists of two and three-body n-n interactions: 

V% c = Mfa - tj) + tjrffo - r,)% - r fc ). (2) 

This interaction @ is supplemented with a long-range 
Yukawa and an effective charge Coulomb interaction 
parts |37!: 
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Zi, Zj are the effective charge of baryons i and j. In QMD 
model, one neglects the isospin dependence of the interac- 
tion. All nucleons in a nucleus are assigned the effective 
charge Z = j^+j^ [37j. The long-range Yukawa force 
is necessary to improve the surface properties of the in- 
teraction. The parameters jj,,tx,t2,t^ in Eqs ©, ([3]) and 
(HJ are adjusted and fitted so as to achieve the correct 
binding energy and mean square root values of the ra- 
dius of the nucleus [37j ■ The Skyrme- type potential can 
be generalized to density dependent form as [13, IH| : 
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with a and j3 as free parameters. These are fixed by the 
ground state properties of normal nuclear matter. The 
third parameter 7 allows to fix the incompressibility of 
the nuclear matter. The two parameterizations are com- 
monly used namely 'soft' and 'hard' corresponding to dif- 
ferent incompressibilities of the nuclear matter. We shall 
use here a 'soft' equation of state with incompressibility 
/C=200 MeV. 



III. NUMERICAL CALCULATIONS AND 
DISCUSSION 



gies of 400, 600 and 1000 AMeV. Figure 1 displays the 
model calculations of the average nucleon density p avg , 
collision rate dN co u/dt, directed transverse momentum 
and momentum anisotropy ratio (Ri SO ) of the par- 
ticipant matter obtained for the reaction 40 Ca + 40 Ca at 
a 'reduced' impact parameter b/b m ax—0.3. The nuclear 
matter density is calculated as: 
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In our approach, the average nuclear matter density p av9 
is calculated in a volume cell of 2 fm with center located 
at the point where two nuclei initially touch each other in 
cm. frame. As expected, one obtains maximum central 
density at the highest considered energy of 1000 AMeV. 
A similar behavior can be seen for the allowed collision 
rate dN co u j dt. The directed transverse momentum 
is another important quantity describing the transverse 
expansion of the nuclear matter. It is defined as: 
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with p x (i) as x-component of momentum, g/j is the ra- 
pidity of i th particle evaluated as: 



_ 1 E(i)+p,(i)c 
Vl 2 m E(i)-p z {i)c 



(8) 



Here E(i) and p z (i) are the total energy and the longi- 
tudinal momentum of the i th particle, respectively. In- 
terestingly at the start of the reaction (for t < 20/m/c), 
{Px lr } turns negative (See figure [He)). This happens due 
to attractive interaction between two nuclei in the be- 
ginning. After the two nuclei collide, strong repulsion 
takes place in the overlap region. As a result, transverse 
flow turns positive and saturates finally. The time evolu- 
tion of the transverse momentum transfer depicts that it 
starts to have finite value when nuclear matter is still in- 
teracting and density is quite high. Beyond 40 fm/c, the 
transfer of momentum in transverse direction is not much 
altered, indicating the expansion of the nuclear matter 
in transverse direction has ceased. One clearly obtains 
more transverse expansion at 1000 AMeV as compared 
to lower incident energies. It means that nucleons with 
larger relative momenta are directed into the transverse 
direction with higher momenta. It is worth interesting to 
compare the degree of equilibrium attained by the partic- 
ipant zone at these bombarding energies. We display in 
the bottom panel of figure [TJ the momentum anisot ropy 
ratio {Riso} for the participant zone, calculated as [8ll22|: 



(Riso) 
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First of all, we study the reaction observables for the 
semicentral collisions of 40 Ca + 40 Ca at incident ener- 



To select momentum anisotropy due to participant zone, 
we impose the condition {-0.5y beam <V< 0.5y beam ) c . m . 
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FIG. 1. The time evolution of 40 Ca+ 40 Ca reaction at incident 
energies of 400, 600 and 1000 AMev and at 'reduced' impact 
parameter b/b ma x=0-3. Results are shown here for: (a) mean 
central density p av9 /p ; (b) average collision rate dN co u/dt; 
(c) directed transverse momentum (pt"~}; and (d) anisotropy 
ratio (Riso) of the participant zone. 



on the event rapidity. The application of this cut allows 
to exclude spectator like components. The full equilib- 
rium corresponds to {Ri SO } value close to 1. One can 
see that equilibration of the participant matter gets sat- 
urated shortly after decompression. This is also the time 
when spectator components start departing from the hot 
interaction zone and we have the uniform directed trans- 
verse momentum, thereafter. This shows that specta- 
tor and participant components hardly interact after this 
time. However full equilibrium is not attained even at 
1000 AMeV. Interestingly, at all three incident energies 
considered here, the anisotropy ratio (Ri SO ) converge to 
same value. This indicates that equilibration of partici- 
pant matter in semicentral Ca+Ca collisions takes place 
to nearly same extent, irrespective of the incident energy. 
These findings are supported by the fact that in Plas- 
tic Ball experiments, nearly the same baryonic entropy 



is produced [4(| at different bombarding energies. The 
Ri SO value is not much altered unlike the mean directed 
transverse momentum (p* r ). It may be mentioned that 
momentum anisotropy ratio is still less than 1 with ap- 
proximately 25 % anisotropy in momentum distribution. 
This also indicates that surrounding spectator compo- 
nents has strong influence on the evolution of participant 
matter. As a result, one can't have always full stopping 
and thermalization even for the central overlap region. 
As far as peripheral collisions are concerned, the par- 
ticipant component, therefore, retains the transparency 
character in the longitudinal direction, and momentum 
distribution of nucleons is far from that of thermalized 
source. This also shows that nucleons are in quite highly 
non-equilibrium environment that can not be described 
with single fireball model [4lj , or oversimplified hydrody- 
namical prescription (42j. 

In [21 we display the momentum distribution of nucle- 
ons for the single event of Ca+Ca collision obtained at 
100 fm/c. The two dimensional plane in the figure corre- 
sponds to longitudinal momentum p z versus total trans- 
verse momentum pr(— yPx+Py) ■ Our model predic- 
tions clearly show the impact of bombarding energy on 
transverse expansion of the nuclear matter. This leads to 
increase in (p* r ) with incident energy, as mentioned ear- 
lier. One can observe appreciable transparency character 
of the nuclear matter in longitudinal direction. This hap- 
pens due to the fact total number of n-n collisions taking 
place in the medium are not able to keep the pace with 
transverse expansion. Owing to this phenomenon, we 
have less equilibrated matter with (Ri SO ) value slightly 
less than unity. Next, we try to understand the stop- 
ping pattern of the nucleons facing maximal and minimal 
number of collisions by analyzing the variation in their 
rapidities at different time steps. Figure [3] shows the vari- 
ation of normalized rapidity y/y beam of single spectator 
nucleon, facing zero collision (n co u — 0) and participant 
nucleon (n co u = 14) followed for the time span of 100 
fm/c. We have displaced here the symbols correspond- 
ing to same rapidity bin but obtained at different time 
steps along the y-axis, though all these points have same 
probability of occurrence. This is done to make the vari- 
ation in the longitudinal rapidity with time more vivid. 
The numbers in the boxes imply the time interval for 
which the nucleon stays in particular rapidity bin. At 
the start of the reaction, both particles are seen at pro- 
jectile beam rapidity yteam ■ With advent of the collision 
process, participant nucleons (suffering 14 collisions) are 
stopped around mid-rapidity with y/yf, eam =0.25. On the 
other side, the spectator nucleons (as shown in the up- 
per panel) still stay in the higher rapidity regime with 
y/y&eoro = l- Again, one can also notice that participant 
nucleon suffers more fluctuations in its longitudinal ra- 
pidity. The spectator particle, however, remains in the 
projectile's rapidity regime for most of the time. The 
temperature reached in HI reactions at intermediate en- 
ergies can be as high as 70-80 MeV 22]. Nucleons feel- 
ing such hot environment would behave differently than 
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FIG. 2. The plots of transverse momentum pr vs longitudinal momentum p z of the nucleons obtained in a single event of 
Ca+Ca collision at 400 (left), 600 (middle), and 1000 (right) AMeV and with b/b max —0.3. Momentum distribution of nucleons 
has been calculated at 100 fm/c. 



those found in relatively cold spectator zones. It is in- 
teresting to follow the space-time characteristics of these 
nucleons facing maximal (i.e. participants) and mini- 
mal (i.e. spectators) number of collisions. This study 
can be of importance to understand characteristics of 
hot participant matter and hadronic interactions as well. 
For this study, we randomly choose the nucleons from 
target and projectile in 40 Ca (400 AMeV)+ 40 Ca col- 
lision. The left panel of 0] shows the time evolution 
of total momentum P„ of individual spectator nucleons 
(P15, P32, -P47, Peo)- Similarly, the evolution of partici- 
pant nucleons (Pg, P39, P49, P71) facing maximal colli- 
sions is shown in the right panel. 

The particles behaving like spectators aren't expected 
to feel the hot environment even during the violent stage 
of the reaction. Owing to this, one obtains minimum fluc- 
tuations in the particle's momentum during the collision 
process. On the other hand, participant nucleons facing 
maximal number of collisions, feel the hot environment 
as indicated by continuously change in their momenta 
during the collision. The participant nucleons, therefore, 
suffer more number of fluctuations in momentum space. 
We also aim to analyze the behavior of these nucleons 
in coordinate space. Figure [5] displays the trajectories of 
these nucleons in x-z plane followed for the time span 
of 300 fm/c starting from the initial contact between 
the projectile and target nuclei. Since nucleons facing 
minimal number of collisions are knocked-out quite ear- 
lier, they traverse nearly straight line paths in 7^3-space. 
Contrary to this, the participant nucleons (as shown in 



the right panel) face large fluctuations along their tra- 
jectories as expected. These particles remain in the hot 
zone for appreciably long time unlike spectator nucleons 
which arc knocked-out of the reaction zone during the 
early stage of the reaction. 

IV. SUMMARY 

Summarizing, we have presented an exclusive anal- 
ysis of the semicentral Ca+Ca reactions at incident 
energies 400, 600 and 1000 AMeV within quantum 
molecular dynamics approach. Our findings showed 
that incident energy of the projectile strongly influences 
the observables viz. average nucleonic density, allowed 
collision rate as well as the transverse expansion of 
nuclear matter. Further, the transparency effect in the 
nuclear medium also seems to get enhanced with the 
incident energy. Interestingly, the degree of equilibrium 
achieved by the participant zone remains insensitive to 
the range of projectile energies considered in this work. 
The full equilibration is not always there even for the 
participant zone. A detailed analysis of trajectories 
traversed by the nucleons depicted that spectator and 
participant nucleons behave differently in the coordinate 
and momentum space. The nucleons facing maximal 
collisions i.e. participant nucleons suffered extensive 
fluctuations in their momentum and spatial coordinates. 
The spectator nucleons, on other hand, depart from the 
hot reaction zone quite earlier and traverse nearly the 
straight line trajectories. 
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